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Abstract 

In this paper, we model the mechanics of a collagen pair in the connective tissue extracellular 
matrix that exists in abundance throughout animals, including the human body. This connective 
tissue comprises repeated units of two main structures, namely collagens as well as axial, parallel and 
regular anionic glycosaminoglycan between collagens. The collagen fibril can be modeled by Hooke's 
law whereas anionic glycosaminoglycan behaves more like a rubber-band rod and as such can be 
better modeled by the worm-like chain model. While both computer simulations and continuum 
mechanics models have been investigated the behavior of this connective tissue typically, authors 
either assume a simple form of the molecular potential energy or entirely ignore the microscopic 
structure of the connective tissue. Here, we apply basic physical methodologies and simple applied 
mathematical modeling techniques to describe the collagen pair quantitatively. We find that the 
growth of fibrils is intimately related to the maximum length of the anionic glycosaminoglycan 
and the relative displacement of two adjacent fibrils, which in return is closely related to the 
effectiveness of anionic glycosaminoglycan in transmitting forces between fibrils. These reveal the 
importance of the anionic glycosaminoglycan in maintaining the structural shape of the connective 
tissue extracellular matrix and eventually the shape modulus of human tissues. We also find that 
some macroscopic properties, like the maximum molecular energy and the breaking fraction of the 
collagen, are also related to the microscopic characteristics of the anionic glycosaminoglycan. 



1 Introduction 

Structural polymers help to maintain the human body shape from large external tractions. Collagen 
occurs extensively in all animals and is the defining structural polymer, existing in the connective 
tissue extracellular matrix (CTs), for example, skin, cartilage and bone. It acts to support life 
against vigorous daily activities. CTs are bridged and bonded by anionic glycosaminoglycan GAGs, 
such as parallel rows of decoran, which are the only molecules in CTs apart from protein fibres that 
can be visualized by an electron microscope [H, Hi 0] ■ A pair of segments on neighboring collagen 
fibrils, which are linked by GAGs chains, is termed a shape module and must deform reversibly 
to preserve the general structure of the organism against various external stresses Q. The axial 
tension is transmitted and opposed by protein fibres while compression is resisted by water-soluble 
polysaccharide GAGs, e.g. chondroitin sulphate. These GAGs transmit forces from the local area 
of molecules to global fibrils by converting compression into disseminated tensile stress The 
mechanics involving the collagen and the elastin fibril is well understood, i.e. they behave mainly 
like a Hookean spring in the low stress limit, but the elasticity of GAGs' molecules is still under 
intensive investigation Q. 
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Numerous computational and continuum mechanical methodologies have been utilized to study 
the mechanics of such CTs 0, HI, S EH 11, 12, 13, 14 1. However, most tend to ignore the microscopic 
details of CTs or assume a simple form of the molecular potential energy in order to make their 
models tractable. In particular, for a continuum mechanical approach, generally a homogeneous 
structure of CTs is assumed. Since the CTs comprise of collagen pairs, here we adopt basic physical 
concepts and simple mathematical techniques to investigate the mechanical properties of a single 
collagen pair. While utilizing a quadratic energy form for the collagen, we incorporate the statistical 
nature of GAGs into our model by utilizing the worm- like chain model, which has both theoretically 
and experimentally been proved to be applicable to wide ranges of bio-macromolecules including 
unstructured DNA, RNA, GAGs and polysaccharide @, EH El, E3, El]- We note that our work 
may form a theoretical basis for experimentalists for their work on collagen pairs. 

GAGs can go through conformational transitions, which are defined by the sudden elongation of 
bio-molecules due to the change in their atomic allocations without an increase in external forces. 
GAGs, e.g. pectinate, have two distinct chair structures, namely 4 Ci and 1 C , 4 [l||, separated by 
an energy barrier of approximately 11 kcal/mol (20| . Moreover, there exists a boat conformation 
1 > A B with an energy level of approximately 5 — 8 kcal/mol above the most stable A C\ chair energy 
2l|. In the low energy configuration, nature prefers the A C\ chair state over other possible states 



because it possesses the minimum energy configuration in comparison to other possible states. 
However, under an external stress, these polysaccharides can undergo two conformational transitions 
reversibly beyond some critical stresses [221 ] . For example, amylose undergoes its first conformational 
transition when the applied force is around 200 pN and its second transition when the applied force 
is around 500 pN. However, some polysaccharides like pigskin DS only go through one transition, 
while other polysaccharides like poly-anionic HA, neutral methyl cellulose and poly-cationic chitosan 
undergo no transitions at all. The reason for these conformational transitions is closely linked to 
the total number of axial or equatorial linkages existing in each pyranose ring (monosaccharide). It 
is evident that the glycosidic linkages may act as levers to generate a sufficient torque to undertake 
the work done, which is necessary to perform ring conformational transitions. Hence only glycosidic 
linkages with axial linkages can generate enough torque to flip the levers beyond a given critical 
stress while equatorial linkages can not. For a more extensive treatment of the conformational 



transitions see Marszalek et al. |22j for details. This phenomena provides a crucial step for a 



molecule to transform from the entropic region into the Hookean regime. Although the above 
interpretation of the kinking conformational changes is largely accepted by most researchers, there 
exists some other interpretations to explain these conformational transitions 23 1. 

This paper is divided into four sections. In section we derive a mathematical model for CTs, 
while in section numerical results and some extension work on CTs are developed. In the last 
section, we present some conclusions. 



2 Theory 

In this section, we consider some simple applied mathematical models for a collagen pair. Since 
CTs contain repeated units of collagen pairs, a collagen pair is denoted by a single unit of such 
repeated collagen pairs, which is illustrated in Figure[TJ While we model collagens utilizing Hooke's 
law, GAGs are modeled utilizing the worm-like chain model, which has carefully taken the entropic 
nature of the molecular chain into account. Due to the symmetry of CTs, we consider a pair of 
collagens with axial and regular GAGs in between, especially bone. In reality, to assemble such 
segments into a whole CTs is a challenging task owing to the complicated molecular interactions 
between fibrils. In addition, to simplify our study of external applied forces, which consist of 
both tensile and compressive forces, only the tensile stress is examined because GAGs can convert 
compression stress into tensile stress and hence the total tensile stress is assumed to be the vectorial 
sum of both tensile and compression forces. 

A well constructed and stable collagen pair is assumed to maintain its structure by the at- 
tachment of GAGs, subject to at least a small perturbation due to molecular interactions between 
collagen-collagen, collagens- GAGs, thermal fluctuations, sudden shocks etc. Suppose that we ap- 
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Figure 1: Collagen pair before stretching, where L\ and L2 denote the collagens 1 and 2 respectively, 
vi . . .vn and u\ ... u at denote the position coordinates of oligomers in L\ and L2 respectively, Ic and Ig 
are the natural lengths of the collagens and GAGs respectively and N is the total number of oligomers 
that could exist in each collagen. 

ply a tensile stress on one end of the collagen 1 (see Figure [2]), where this perturbation alters the 
mechanical structure of the collagen pair. 

The displacements between Vi + i and Vi are denoted by for i = 1, . . . , N — 1 where N is the 
total number of oligomers, which are defined by a segment of the collagen pair with two GAGs 
attached at its ends. Likewise for the collagen 2, the displacements of iti+i and Ui are denoted by 
5j. Schematic diagrams of an un-stretched and a stretched collagen pair are shown in Figure [T] and 
Figure [2] respectively. We postulate the potential energy of collagens utilizing Hooke's law, which 
reads 

V = \us\ 

where k is a spring constant and 8 is the extension. Given the potential energy form, it is straight- 
forward to show that the potential energy of a collagen pair, V c , is given by 

i=i v 7 i=i v ' 

where N, k c and £ c are the total number of oligomers that exist in the collagen, the spring constant 
of the collagen, and the natural length of a collagen respectively. Moreover, Si and Aj are defined 
by Si = Uj+i — Ui and A, = w^+i — Vi for i — 1, 2, . . . , N — 1 respectively. In addition, we postulate 
the potential energy of GAGs utilizing the worm-like chain model. It is well-known that the 
interpolation force-extension formula for the worm-like chain model 24] is given by, 

'-^{i+^-r-i}. 

where /, ks, T, A, z, L denote the applied force, Boltzmann's constant, the absolute temperature, 
the persistence length, extension and the total contour length of the GAG respectively. Then, the 
potential energy of all GAGs, V g , can then be obtained by integrating Eq. [2] with respect to z and 
summing up the total number of GAGs, N, to yield 
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Figure 2: Collagen pair after stretching, by an applied force F assumed to be acting on L\ causing 
an induced force F' in L2, while s denotes the offset length between L\ and L2, Ai and b\ denote the 
length of the first oligomer in L\ and L2 respectively and so on. 




where k g and £ g denote (fcsT) / (AL) and the natural length of GAGs respectively. We can then 
relate it, and Vi to Si, Aj and s, which is the off-set between fibrils (See Figure [5]). Notice that s is 
a function of / as the collagen pair starts to slide away with respect to each other, subject to the 
external force. Hence, Ui and vi are geometrically related by 

Vl = iti + s, 

v 2 =u 2 + s + (Ai - £ C /N) - (<Ji - t c /N) = u 2 + s+ (Ai - 81), 



Vi = m + s + ^(A fc - 4), 
fc=i 

N 

VN = u N + s + ^2(A k - 8 k ). (4) 
fc=i 

After we relate the kinematics between the collagen pair, we can simplify Eq. [3] in terms of the 
displacements <5j and Aj, giving 
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(5) 



Since the elongation of collagen 1 causes the elongation of collagen 2, we can relate Si and Aj by 
the relative displacements ej. That is 



Aj = «i + ei. (6) 

Assuming statistical equilibrium, we let Si = S, A» = A and ej = e for all i. Hence, Eqs. [T] and \5\ 
reduce to 




4 1 - L 



k g \J{s + ief + {i g f-i g \ \. (7) 



Since each GAG has its own maximum length K and the displacement between un and corre- 
sponds the maximum length of GAG given in the collagen pair, which must be smaller or equal to 
K. Given that, we have 



\v N - u N \ = yj(s + NeY + {l g y < K, (8) 

from which, we show that the maximum number of oligomers N max allowed in a given collagen 
satisfies the following inequality 



e e 

where |_ J denotes minimum integer value of the enclosed real number, e.g. [4-3J = 4 respectively. 
Note that Eq. [9] ignores the maximum strain that can be held by a collagen, and while e is a rather 
abstract quantity, it can be expressed in terms of s by minimizing the maximum potential energy 
E m ax of the collagen pair with respect to e at equilibrium (see Appendix for details). At first 
glance, s is insignificant under small natural external forces, i.e. thermal fluctuations and internal 
molecular interactions. However, s becomes crucial when we consider the collagen pair under large 
external tractions. One interesting thing about the above equation is that it limits the possible 
number of oligomers that can exist for a given stable collagen and it is model independent because 
it arises solely from a geometric point of view. Also, the inequality reveals the importance of GAGs 
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in relation to the growth of the collagen pair. Hence, the longer the maximum length of GAGs, the 
longer the structural collagen pair can be. Although the inequality does not prove the possibility of 
the existence of ring conformational transitions of GAGs, it does provide evidence that the existence 
of conformational transitions increases K and hence encourages the growth of the collagen pair. 
Further, the higher the effectiveness of the force transmission of GAGs between fibril and fibril, the 
lower the value of e, and hence longer the collagen can be. 

The induced force, F', in collagen 2, generated by collagen 1 and the relative displacement, e, 
of the system, can be determined to be 

F' = F-k c e, e=^(Li-L 2 ), (10) 

where L\ and Li are the total molecular contour length of the collagen 1 and 2 respectively. F' 
and e can hence be determined experimentally by knowing k c , F, N, L\ and L^. In particular, e 
can be measured to find N max . Given N max , we can obtain the maximum molecular energy for a 
collagen pair, E max , namely 




The maximum molecular energy E max grows quadratically with the extension 5 and is related 
linearly to N max , which is given by Eq. and determined by K and e. Hence, the existence of 
flexible GAGs with long maximum length K increases the number of oligomers N max , which in 
return tightens the system dramatically. In addition, assuming a large N rnax and a small e, the 
modulus of a collagen pair, fc t , can be determined easily by the second derivative of E max with 
respect to S. Upon performing two differentiations, we approximate the modulus of the collagen 
pair, A;*, as 

kf ~ k c N max . (12) 

Again, k t depends linearly on N max , where the rest of arguments are very similar to the discussion 
for the maximum molecular energy above. 

3 Numerical results and analysis 

In this section, we carry out a numerical analysis on the results derived in the previous section, and 
finally we examine the breaking point of a collagen pair. Firstly, given Eq. 111! we assume the value 
of the parameters given in [7j, namely k c = 0.2 GPa, k g = 0.02 GPa, l c = 100 /an, £ g — 0.02 nm, 
L = 0.2 nm and 3 maximum numbers of oligomer are examined, namely N max \ — 1 K, N maX 2 = 10 
K and N max ^ = 0.1 M. Further, without loss of generality, the offset s is assumed to be zero, 
which corresponds to the tight molecular interaction between fibrils or the collagen pair under a 
small traction. The potential energy of GAGs is neglected due to its insignificance in comparison 
to the potential energy of the collagen pair that is demonstrated in Figure [31 where e has been 
linearized with the extension 6. The numerical result shows that the potential energy of GAGs 
utilizing the above parameters is 6 orders of magnitude smaller than the potential energies of the 
paired collagen for all three cases we considered. Given that, the potential energies of the collagen 
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Figure 3: Potential energy of GAG, V g , versus extension, 5 ranging from to 100 nm, for N max \ = 1 



K, N max2 = 10 K and AT, 



max3 



0.1 M respectively. 



pair versus the extension S, ranging from to 100 nm for N max i, N maX 2 and N maX 3, are plotted 
together in Figure |U 

Note that the maximum molecular energy increases quadratically when the extension increases 
linearly. In addition, i5 is determined by the applied force F and the relative displacement e. For 
example, for a constant A, 8 achieves its maximum value when e = 0. Further, E max increases 
linearly with N max , which in return depends positively on K but negatively on e. The toughness of 
a collagen pair, k tl is defined in the previous section as the second derivative of E max with respect 
to S, and hence the higher the slope of E max with respect to the extension 6, the stiffcr the collagen 
pair is. In conclusion, given a collagen pair, it is tougher whenever GAGs have a higher value 
of K and a higher ability in transmitting forces between fibrils. However, the existence of GAGs 
contributes almost nothing to the toughness of the collagen pair but is significant to maintain the 
stable structure and boost up the toughness of the collagen pair, which are consistent to the results 
that an increased PYD or DPD ratio (the most abundant mature GAGs in bone collagen) is related 
to the increased compressive strength in bone (CTs) [1^, [26|, [U [H, [2§| but has no huge effect on 
toughness or ductility of bone [3(J EH, H3, [33| . 

We end this section by extending concepts that we developed above to determine the point at 
which a collagen pair will break. For the study of a collagen pair, an important parameter is the 
fraction of the collagen pair, XRi m which the breaking point occurs (see Figure [5]). The breaking 
point is assumed to occur when the bridge between u\ and V\ is finally broken (see Figures Q] and 
[5]). From Figure we write 



t max 
L 1 



XR = sr + 



(13) 



where L™ ax and L™ ax 



are the critical lengths of collagens 1 and 2 respectively and sr is the critical 
off-set of the system. After some re-arrangement, we find that 



XR = LT X - s R . (14) 

In addition, we know that when the collagen pair is about to be torn apart, N max — 0, from Eq. 
we deduce 



sr. 



{Kf 



(15) 
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Figure 5: The breaking fraction of the collagen pair 
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Therefore, upon comparing Eqs. 1141 and 1151 we have 



X r = LT X ^{Kf - (l g y. (16) 

Notice that the above equation links the macro-quantity \R to micro-quantities such as L™ ax , K 
and i g . Once again, the larger the possible length of GAGs, the smaller the breaking fraction of 
the collagen pair, which implies that we need to put more forces in the system to tear the collagen 
pair apart, which is consistent to the theoretical result obtained by [HI]. 



4 Conclusion 

In this paper, we utilize simple applied mathematical modeling techniques, for example Hooke's 
law, to describe the structure of a collagen pair. We find that for maintaining the stability of a 
collagen pair, the maximum number of oligomers that can exist in a given collagen depends on 
the maximum length of GAGs and the effectiveness of the GAGs in transferring forces between 
fibrils. That is, a collagen can grow longer for the longer GAGs and the higher effectiveness of 
GAGs in transferring forces between fibrils. This concept can then be extended naturally to the 
toughness and the breaking point of a collagen pair, which are found to be intimately related to the 
characteristics of GAGs. In addition, the possible conformational transitions of GAGs strengthen 
the whole structure of CTs and as such CTs maintain our body shape from any large external 
traction. We note that our work does not give the complete picture where a more sophisticated 
theoretical investigation of the interactions between collagens and GAGs needs to be undertaken 
and extend our model to investigate the mechanical properties of CTs. In addition, we sincerely 
refer our readers, who are interested in utilizing computing simulations and continuum mechanics 
to model the macroscopic structure of CTs, to those references listed in the introduction. 

Appendix 

From Eq. Ill) if N max is sufficiently large, then we can approximate the summation by integration. 
Upon utilizing equation Eq. [9j we obtain 



Em.ax — ~ k c 



1, /$-S 



<5- 



$ — s 



4>-s , 



d.r 




(s + exf + £ 2 g -£ g 



1 - 
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-1 2 
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-k a {^j{s + e X y+P g -t g }dx, 



(17) 



where $ = \ K 2 — I 2 . If we substitute s + ex — £ g tan0 and perform the integration, we obtain 



1 /$-s 

^max — 2 \ 



$ - S 



$ - s 



H(s) 



(18) 



where H{s) = {l/2)k g £ g a(s) + (L 2 / A)k g t g b{s) - {L / A)k g i 2 g c(s) and 
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a(s) = 2 tan 9 + (1/3) tan 3 9 - tan 9 sec 6» 

-Msec^ + tan^Cr!!/^. 

L In [tan | + l] L In [tan § - l] 
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(19) 



We note that the energy form given in Eq. [18] might be utilized to carry out computer simulations 
or mathematical modelings for CTs in a more accurate way. To find the minimum value of E max 
with respect to e, we require dE rnax /de = 0. Given that, we have 



\ 



S 2 



H 



fc c (*-s) 



<t> — s 



1 - 



<E>— s 



(20) 



Upon knowing s, e can be easily solved from the above equation. 
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